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The  work  discussed  in  this  report  was  performed  under  Work  Unit  72330416,  Self-Organizing  Control/RPV 
Augmentation,  This  work  suDPlementaAMRL  Technical  Report  AMRL-TR-73-130>  TERCOM  Performance: 
Analysis  and  SimuIationT by  presenting  the  effect  of  two  common  errors,  heading  at»d  ground  speed,  on  the 
TERCOM  model  derived  in  that  report.  Improvements  of  terrain  generation  have  atso  been  incorporated  so 
that  realistic  non-Gaussian  terrain  may  be  treated.  The  analysis  and  simulation  work  was  performed  in  the 
Mathematics  and  Analysis  Branch  of  the  Biodynamics  and  Bioengineering  Division,  Aerospace  Medical 
Research  Laboratory  during  the  period  March  1975  to  August  1977. 
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INTRODUCTION 
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Thf  terrain  contour  matching  tTKRC'OM'  uaviKotion  ayatem  ia  propoainl  aa  an  autoinatir  updatiOK  lievut*  to  Ih> 
uaetl  aa  a aupplement  to  atandard  inertial  navigation  ayatema.  The  aircrart  uainu  TKRlHtM  aainplea  the 
terrain  amplitudea  beneath  ita  flittht  path  by  meana  of  a radar  altimeter  and  atorea  a atrinn  of  N aarnplea  in  a 
temporary  on-lHU»nl  computer  memory.  Thia  vector.  itm\poa»>d  of  the  N-terrain  aarnplea.  ia  then  compannl  with 
all  N wmponent  amplitude  vectora  extracted  t'rom  an  amplitude  array  in  permanent  memory,  repreaentinjs 
the  topography  of  the  terrain  over  which  the  aircroR  la  expect  eti  to  Im  flying  The  Imat  match  between  aamplinl 
and  atoreil  vector  determinea  the  kihuiiuI  track  over  which  the  aircraft  haa  juat  flown  and  hence  ita 
KeoKraphical  poaition  Thia  ayatem  i-eipiiiva  that  the  inertial  ayatem  Ih'  accurate  enouKh  to  place  the  aircraft 
inaide  the  Imundanea  of  the  terrain  map  atorinl  permanently  in  the  on-hoanl  rtunputer  \ typical  mi.aaion 
would  ctmaiat  of  a number  of  inert  ial  leija  with  a TKKlHlM  poaition  check  at  the  end  ol  each  to  aaaure  more 
accurate  navigation  than  wuld  be  achievtHl  with  the  inertial  ayatem  alone  I'he  alnonthina  pi-opoawl  for  uae  in 
matching  the  atortnl  terrain  vtn'tor  with  the  sampled  terrain  vector  are  the  M.\n  uninimnm  average  distance' 
cla.saifier  and  the  MSP  iminimuni  mean  aijiiartHl  distance'  claasifter  The  M.M'  claaaifter.  which  is  preterred 
biwauae  of  itaiMinputational  simplicity . was  studied  in  a pivvions  report  iref  1'  and  was  shown  to  j;ive 
essentially  the  same  performance  as  the  MSP  classifier  under  all  conditions  studied  under  ret  I 


TERRAIN  SAMPLING  ERRORS 


In  the  previous  report  (ref  I)  the  terrain  amplitude  sampled  by  the  aircraft  was  assumed  to  have  had  the  same 
((eographical  coordinates  as  some  subset  of  the  terrain  amplitudes  stored  in  the  on-board  computer  memory. 
Such  a condition  is  shown  in  Figure  1,  in  which  a sample  terrain  section  is  shown.  The  h axis  designates  the 
terrain  amplitude  above  (or  below)  some  reference  plane.  The  amplitudes  representing  this  terrain  that  are 
stored  in  the  on-board  computer  are  denoted  by  the  length  of  the  vertical  line  between  the  reference  plane  and 
the  heavy  dots  at  the  intersections  of  the  X and  Y grid  lines.  The  row  of  crosses  running  across  the  middle  of 
the  terrain  represent  the  location  of  the  amplitudes  sampled  by  the  aircraft  as  it  flies  over  the  real  terrain. 
Gaussian  noise  was  added  to  the  stored  terrain  to  account  for  aircraft  sampling  errors,  errors  in  stored  terrain 
altitude  values  and  all  other  possible  errors.  This  was  certainly  an  oversimplification,  since  all  noise  sources 
were  assumed  additive  at  all  levels.  This  assumption  is  not  true  for  a combination  of  heading  and 
ground-speed  errors,  as  shown  later  in  the  report. 


Figure  1.  Terrain  sample  showing  stored  amplitudes  and  sampled  amplitudes  for  original  TERCOM 
analysis. 

A more  realistic  approximation  is  shown  in  Figure  2,  where  the  aircraft  sampling  of  altitudes  may  occur 
anywhere  on  the  x,  y plane.  Figure  2 shows  a possible  sampling  track  in  which  the  distance  between  aircraft 
sample  points  ( AR*)  is  still  equal  to  the  distance  between  the  .sample  points  stored  in  memory  ( AXt'  and  the 
track  runs  parallel  to  the  axis.  Figure  .1  shows  a condition  where  aircraft  ground  speed  is  slightly  faster  than 
the  programmed  speed,  so  the  distance  between  aircraft  sample  points  (AR,»)  is  greater  than  .\Xt.  If  ground 
speed  were  lower  than  the  programmed  .speed,  AR*  would  be  le.ss  than  AXf  Figure  4 shows  an  example  of  a 
course  heading  error  <r.  The  track  is  no  longer  parallel  to  the  X axis,  but  AR*  “ AX,,  so  the  ground  speinl  is 
correct.  A final  condition  where  both  AR*  and  it  vary  was  studied  so  that  the  interaction  betwinm  the  two  could 
be  investigated. 
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Figure  2.  Terrain  sampte  showing  stored  amplitudes  and  sampled  amplitudes  with  no  groundspeed 
or  heading  error. 


Figure  3.  Terrain  sample  showing  stored  amplitudes  and  sampled  amplitudes  with  a groundspeed 
error  but  no  heading  error. 


TERRAIN  GENERATION 


Experience  with  terrain  generation  gained  in  previous  work  showed  that  typical  real  terrain  arrays  contain 
both  flat  and  hilly  areas  tref  11,  Amplitude  histograms  computed  from  a large  area  of  real  terrain  appeared  to 
be  nearly  Gaussian.  However,  amplitude  histograms  computed  from  smaller  ureas  of  this  same  terrain  array 
were  non-Gaussian.  (The  distance  between  sample  points  on  these  arrays  was  400  ft.)  The  non-Gaussian 
behavior  of  the  smaller  terrain  samples  and  visual  inspection  of  these  same  samples  indicated  that  more 
realistic  terrain  could  be  generattnl  from  a non-stat ionary  amplitude  distribution.  The  comparison  of 
simulation  runs  for  real  and  artificially  generated  terrain  in  ref  1 showed  that  TERCOM  performance  over 
artificial  terrain  was  the  same  as  for  real  terrain  if  a certain  parameter  ratio.  (T/UTy,  was  the  same  for  both 
types  of  terrain.  The  parameter  it,  is  defined  as  (r,*  E{(X,  - X,,,)*}  where  E denotes  expectation  and  the  X, 

are  the  terrain  sample  points,  ir,-  is  just  the  variance  of  the  sample  amplitudes,  it,  gives  a measure  of  terrain 
roughness,  but  it  was  found  that  <T,;iTy  was  more  informative  than  ij,  alone.  A small  n,  (ri-  indicates  smooth 
terrain  with  little  variation  between  the  sampled  amplitudes,  but  possible  large,  slow  fluctuations  in 
amplitude  over  the  entire  area.  A large  (cr,  (rr>  indicates  large  variation  between  at^acent  sample  points  but 
possible  small  variation  in  overall  amplitude.  It  was  found  that  iT,UTy  values  in  the  range  required  could  not  b<> 
produced  by  artificial  terrain  generated  by  a stationary  Gaussian  process.  The  first  methiHl  used  to  introduce 
nonstationarity  was  simply  to  generate  a Gaussian  distribution  of  terrain  amplitudes  and  then  to  set  all 
amplitudes  in  a number  of  adjacent  rows  of  the  array  equal  to  zero  to  produce  fiat  areas.  The  number  of  rows 
set  equal  to  zero  was  acfjusted  until  the  value  of  iT,UTy  matched  that  of  a particular  set  of  real  terrain 
amplitudes.  A different  approach  was  clearly  required  to  remove  the  amplitude  discontinuities  present  in  this 
procedure.  The  algorithm  finally  derived  involves  the  generation  of  a (Jaussian,  correlated-amplitude  array 
from  the  method  developed  by  Moshman  (ref  31  and  described  in  Appendix  B of  ref  1.  and  the  multiplication  of 
this  array  by  a gain  function  that  varies  with  distance  in  one  direction  across  the  array. 

The  gain  function  is 

GiXl=  1 

(1  f X X,v^ 

In  terms  of  terrain  sample  points,  X,.  where  i = 1.  (>4. 

(;,X|)=  1 where  X„  = 32. 

|1  ( (X,-irX.,p 

The  value  of  er/  irp  is  adjusted  by  changing  X„.  the  coordinate  where  G(X,1  = 

This  IS  the  form  of  (i(  XI  used  by  the  computer  program  in  appendix  A.  A typical  terrain  array  generated  by 
this  program  is  shown  in  Figure  5,  An  amplitude  histogram  of  this  array,  shown  in  Figure  ti,  displays  a 
non-Gaussian  shape.  However,  ten  arrays  were  generated,  and  when  all  were  normalized  to  have  the  same 
mean  amplitude  of  zero,  the  histogram  resulting  from  the  combination  of  all  ten  is  essentially  (uiussian,  as 
shown  in  Figure  7. 


where  X„is  the  X coordinate 
of  the  center  point  of  the  array. 
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Ki^jurt'  5.  Typical  terrain  array. 
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KiKure  M.  .Amplitude  hi.ttoKrani  of  terrain  in 
Kijiure  5. 


Fij(ure  7.  .Amplitude  histoitram  v»f  all  10  terrains 
combined. 
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FLIGHT  SIMULATION  PROGRAM 


This  program  was  soparato  iVom  the  terrain  jjeneration  program  aiul  used  ten  pre-prepared  terrain  arrays  to 
simulate  TKRl'OM  runs  with  variations  in  ground  spivd  and  course  heading  The  program  read  in  a 
particular  t'4xd4  terrain  array  and  set  up  two  repivsentations  of  the  array  in  computer  memory  One 
reprvsentation  larray  S'  was  servtni  as  the  array  stored  in  the  aircratVs  on-hoard  memory  Noise  of  a selected 
level  was  addini  to  this  array,  point  hy  point,  to  simulate  ernu's  intnKluced  in  transcrihin>j  terrain  amplitudes 
Irom  small  scale  aerial  photographs  or  maps  iref  1?'  I’his  noise  ma>  he  correlaltnl  or  uncorrelate<l,  hoth  i vpes 
were  testtnl 

I he  second  representation  (array  K'  was  us«>d  Ic  simulate  the  actual  terrain  over  which  the  aircrart  was  tlviiij; 
The  radar  altimeter  values  were  obtained  from  tins  array  I'he  track  len>;th  of  each  rKKl'llM  run  was  4S 
points  out  of  the  availahle  h4  Kollowinj;  the  information  in  ref the  distance  iH't  ween  points  was  set  tHjual  to 
400  rt  Kach  such  distance  will  Iv  refeiTtnl  to  as  a cell  lenfjth  I'his  makes  the  average  4S  cell  track  len^tth 
tH)ual  to  Ith'.’t'O  It  or  alnuit  0 0 miles  IVack  starting  coordinates  were  chosen  at  random  hut  limittHl  to  assure 
that  the  entire  4S  point  track  would  lie  w ithin  the  ti4x04  array 

The  runs  proceeded  in  a dinvtion  parallel  to  the  \ axis  It  was  assunu'd  that  terrain  sampling:  noi.se  in  the 
radar  altimeter  .samples  from  array  K would  Iv  much  smaller  than  emus  present  in  the  on  hoard  computer 
representation,  array  Siref  2'  Therefore,  no  noise  was  addtHi  to  the  track  ohtaiiu'd  from  array  K 

.•Vs  mentioned  earlier,  this  simulation  allows  the  track  to  hefjin  at  any  X and  Y co«>ixliiiale  ithat  would  allow 
the  track  to  Iv  completely  within  the  array'  and  not  just  at  the  discrete  points  represented  in  the  on-hoard 
terrain  array  Moth  terrain  arrays  represented  in  the  Hijilit  simulation  pro>;ram  exist  only  asdiscrete 
amplitudes  at  coordinates  x,  and  y,,  which  I will  call  1 and  d,  n'sp«»ctively  The  metluHl  u.sed  to  extrapolate  the 
terrain  lH'tw«vn  the  I and  .1  c(HUxlinates  in  array  K.  for  use  as  the  radar  altimeter  sampled  values,  is  verv 
simple  The  three  sets  of  d.  .1'  ivordinate  v>airs  nearest  to  each  (\,  Y'  rtHv.dinate  pair  are  determined  hy  the 
pntgram.  and  a plane  is  defined  hy  the  tlm>e  sets  of  terrain  amplitudes  .\(l,  tl'.  The  terrain  amplitude  at  i\,  Y' 
calKHl  Hi\.  V ' IS  then  defined  as  the  distaiux'  hetween  this  plane  and  the  H I'  plane  .-Vn  illustration  of  this 
linear  interpi'lation  met  tun!  is  shown  in  Ki>;iux'  S In  a typical  track  >;enerat«Hl  hy  the  pro>;rani.  none  of  the  X.Y 
ciHirdinate  pairs  coincide  with  an  I..I  cooixlinate  pair 


All  ,j»n 


Figure  H.  Terrain  amplitude  extrapolation  method. 
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1\)  siimilau'  Kioutul  «‘rrors,  a paramotcr  ratio  V V„  was  oiUorod  into  tlio  pro»;ram  \'  was  tlio  ti  lu-  uroumi 

sjuH'd  while  V,,  was  tho  pn'jjramnuHi  01  expiYtinl  ^troiiiui  speed  at  whtcIwXK  .\\,  1 lhema\imum\  \„i'atio 

iistxi  was  t 1 and  the  miiiiiiuiiii  was  i*.  represent iti>;  ahovit  a ♦ U>‘'!  change  in  (jroiind  speed  1 his  prinlneed  a 
ma\iimim  traek  len):th  ot'4S  \ 1 1 = 5‘J  S cells  This  lengthening;  of  the  sampled  track  constrained  the  \ 
ciHUxl  mates  of  the  randomly  selected  st  art  in>;  points  to  X • 1 1 in  order  that  t he  ent  ire  IS-pi'int  1 1 ai  k should  I it 
within  the  array 

I'he  final  elYect  w as  the  introtluction  of  course  heading;  errors  I’he  maximum  an>;ular  error  considered  w as  S 
Phis  further  limited  the  random  startiiii;  ciH'rdiiiates  alon>;  the  y axis  It  the  maximum  possible  track  length 
was  .“Sd  cells  and  the  maximum  angle  was  the  projtH-tioii  of  this  track  on  the  y axis  was  od  sin  S , which  is 
abn'iit  7 4 cells  I'he  y axis  starting  ciHU'dinates  were  thus  limittHl  to  the  region  d • y - oiv 
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SIMULATION  RESULTS 


TERCOM  PERFORMANCE  WITH  SYSTEM  NOISE  ONLY 

The  terrain  generation  and  track  placement  methods  were  changed  sufficiently  between  this  report  and  the 
previous  one  (ref  1)  that  it  is  worthwhile  to  compare  some  results  from  ref  1 to  the  new  data.  Most  of  the 
performance  results  in  the  previous  report  were  plotted  in  terms  of  the  number  of  correct  position  fixes  vs. 
iT/JfTf,  (o-N  is  the  standard  derivation  of  the  additive  noise).  Since  none  of  the  tracks  in  the  present  simulation 
fall  on  the  exact  I and  J coordinates  of  the  stored  terrain,  a correct  position  identification  will  be  defined  as  a 
fix  (in  terms  of  I and  J)  within  one  cell  length  of  the  track  starting  coordinates.  The  values  of  tT/Ja-f  for  the  ten 
terrains  generated  range  from  . 1 to  .33,  so  the  results  were  compared  to  a performance  curve  from  ref  1 which 
had  been  computed  on  terrain  with  a <T,/iTr  of  .2.  This  comparison  is  shown  in  Figure  9 for  the  case  when  noise 
added  to  the  terrain  was  uncorrelated.  The  solid  curve  represents  the  data  of  ref  1 and  the  points  represent  the 
new  data.  Each  point  represents  the  average  of  ten  runs  on  a particular  terrain  array  at  the  level  of  tT/JtTs 
indicated.  Clearly  there  is  almost  no  change  between  new  and  old  data.  Starting  the  track  at  points  between 
the  1,J  geographical  coordinates  has  no  effect  as  long  as  the  ground  speed  and  course  heading  are  the 
programmed  values.  Another  method  of  looking  at  the  data  is  shown  in  Figure  10  where  the  average  miss 
distance  is  plotted  as  a function  of 


F'igure  9.  Percent  correct  responses  (within  I cell  of  correct  location)  for  tn  = 0 tt  = 1024  cells. 


It  is  evident  that  for  low  <t,  ityf  values  (high  noise  Itwels)  the  error  in  location  can  be  as  large  as  (iOOO  to  SOOO  ft 
However,  for  values  of  (r, /<rN  greater  than  1 .fj  we  have  an  almost  perfect  correct  fix  score.  A similar  plot  is 
shown  in  Figure  1 1,  where  the  noise  adde<i  to  the  terrain  is  correlated  and  has  the  same  correlation  length  as 
the  terrain.  In  this  case,  the  miss  distance  does  not  descend  as  rapidly  as  it  does  in  Figure  10.  and  perfect  fix 
capability  is  not  reached  until  ir//<r,v  - 2.5.  The  average  miss  distances  at  lower  <r//<r^  values  are  also  higher 
than  tho.se  shown  in  Figure  10. 


TERCOM  PERFORMANCE  WITH  VARIATION  IN  GROUND  SPEED 
AND  COURSE  HEADING 

Th»>  ortwt  of  Kmund  8}u*o<l  varialion  waa  carrioci  out  iintior  conditiouH  whoro  «r,.(rN  was  than  100.  aiul 

was  jfroatt'f  than  10,  thorohy  ininiini/.inK  noiao  otTwta  IVoin  aourooa  othor  than  urouiul-apowt  variation. 
This  aiinulation  aaatimoa  that  th**  airoraft  ia  aainpliiiK  terrain  ainplitiuioa  at  a fixinl  rate  ao  that  variationa  in 
Kround  apewl  about  V,,  will  rauae  the  Kt'^Kfaphieal  diatance  ludween  aainple  pointa  to  vary  about  the  expected 
4t)t)-ft  intervala.  The  effect  of  a constant  j{«'ound-apeed  error  ia  shown  in  Finure  1'.!  where  an  asymmetry  in 
performance  la  evident  about  the  value  V V,.  « 1 It  ia  seen  in  Kinure  12  that  the  averanc  and  the  variance  of 
the  miaa  distance  ({row  more  rapidly  for  \'  \\  - I than  for  \'  V„  ■ I an  unexpwted  result  Fi({ure  Id  shows 
the  increase  in  averajte  miss  distance  as  the  course  headin((  error  increases  from  0"  to  H".  A({ain,  each  point 
represents  10  trials  on  a particular  terrain,  and  the  mean  is  computed  from  100  (rials  distributed  over  10 
terrains  Mere,  ir,  is  ({r»>ater  than  2;  larne  enou({h  to  ensure  an  averaite  miss  distance  under  1 cell  len(;th  for 
(V  heading  error. 


TfHRAIN  CORRtLATiON  LENGTH 
TL  • 1024  CELLS 

NOISE  CORRELATION  *0 


90  M 10  I 05  II 

V/Vo 


Fiifurt'  12.  Variation  in  averaite 


miss  diatanre  with  ain'raft 


N|M>ed  for  no  hoadiiiK  err«>r. 


II 


TRACK  ANGLE  IN  DEGREES 


Kiguro  13.  Av»»r«j{o  miss  distunrc  ns  ii  function  of  truck  iinttlc  for  • 2. 


Finally,  Fi^jiircs  1 •!  and  l.S  show  the  cITcct  orKround-speed  errors  lor  3 dillerent  course  headnifi  errors  h ij;ure 
It  shows  the  results  fora  track  len>{lh  of  4H  cells  while  Ki>;ure  Ih  shows  performance  with  a shorter  track 
length  of  113  cells. 

The  increased  miss  distance  for  ground  speeds  lower  than  V„  is  more  evident  for  lar^;er  course  lieadm^  errors 
anil  smaller  track  lengths. 
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CONCLUSIONS 


This  siimilation  provided  a more  realistic  simulation  ol  TKKt'OM  than  the  method  used  in  the  previous  report 
irel  1 1,  The  simulation  showed  that  even  when  the  radar  altimeter  samplint;  points  fell  between  the  X and  Y 
coordinates  of  the  amplitudes  stored  in  the  on-hoard  computer,  the  per  cent  of  positions  identified  correctly 
remained  about  the  same  function  of(r/  ox  as  those  of  ref  1 

Apparently,  there  is  sutricient  correlation  between  sample  points  in  the  terrain  model  so  that  correct  fixes  can 
he  made  with  probability  • .ft,')  for  it,  15  if  the  additive  system  noise  is  uncorrelated.  With  correlated 
noise  having;  the  same  correlation  len^'th  as  the  terrain,  we  saw  that  the  95‘(  correct  fix  ran>'e  is  it,  it^  ■ 2.0. 
This  result  is  in  agreement  with  the  findinj;s  of  ref  2 that  correlated  noi.se  caii.ses  a hij;her  error  rate  than 
uncorrelated  noise.  For  lonj'er  track  leiif't  hs,  it  is  expected  that  the  it,  ir;.(  ratio  for  95' i correct  fixes  would  he 
lower. 

TKRl'DM  may  he  used  with  unmanned  vehicles  such  as  lonn-ran^je  Kl’N's.  I’erformance  of  precision  attack  or 
reconnaissance  missions  will  require  remote  operator  intervention  during  the  target  approach  phase  of  the 
tliuht.  Present  control  philosophy  envisions  a video  feedback  to  the  operator.  enahliiiK  him  to  ident  ify  the 
tar^jet  or  specific  landmarks  and  to  make  course  corrections  or  weapons  launches  as  required.  However,  due  to 
the  enemy  electronic  detection  and  the  suh.sequent  jamming  threat,  it  is  felt  that  video  transmission  time 
should  he  held  to  a minimum,  and  consequently  video  frame  rates  may  he  le.ss  than  normal.  .Ml  this  implies 
that  the  KPV  operator  should  find  his  vehicle  close  to  the  target  when  viileo  transnii.ssion  begins,  or  he  may 
not  he  able  to  locate  it  at  all.  Studies  underway  at  AMRl.  RHM  indicate  that  even  with  a narrow  field  of  view, 
say  29'’.  and  an  altitude  between  500  to  1000  ft,  a positioning  error  of  around  1000  0 still  gives  the  RPV 
operator  enough  time  to  maneuver  the  vehicle  over  the  target  or  to  release  weapons. 

If  we  then  consider  SOO  tt,  or  2 cell  lengths,  to  he  a maxinuiin  acceptable  average  position  error,  we  see  from 
Figure  1-1  that  this  error  occurs  for  a course  heading  deviation  of  2"  from  the  expected  cour.se.  Figure  1 1 shows 
that  this  criterion  is  also  exceeded  for  a ground-speed  error  of  +5''; , The  interaction  between  the.se  two  error 
sources  shown  in  Figure  14  tells  us  that  if  a course  heading  error  of  2"  is  allowtHl,  the  ground-speed  error  must 
be  limited  to  ±3'f  if  we  are  to  meet  our  800  ft  criterion.  The  ground-.speed  error  in  the  simulation  occurs 
hwause  TKRCOM  is  assumed  to  sample  altitudes  at  a fix  time  rate,  while  the  actual  aircraft  ground  speed 
may  he  somewhat  higher  or  lower  than  expected.  If  a doppler  radar  sy.stem  can  he  u.sed  to  measure  ground 
speed  with  sulTicient  accuracy,  the  ground-.speed  error  could  he  reduced  by  altering  the  sampling  rate  of 
TP^RrOM.  However,  even  if  ground-.speed  errors  are  eliminated,  the  error  in  the  course  heading  must  he  kept 
to  below  2"  for  the  48-cell  track  length.  While  the  48-cell  track  length  and  the  t>4x04-cell  terrain  matrix  are 
reasonable  for  an  operational  system  (ref  2>,  other  longer  track  lengths  and  larger  array  sizes  may  he  preferred 
for  greater  positioning  accuracy. 

Neverthele.ss,  the  calculations  show  that  the  accuracy  of  the  system  depends  most  heavily  on  the  course 
heading  error.  This  in  turn  is  related  directly  to  the  accuracy  of  the  inertial  navigation  sy.stem,  and  cannot  he 
corrected  by  doppler  radar,  as  can  the  ground-speed  errors.  Thus,  for  any  system  using  only  inertial  and 
TF'RCOM  navigation  (with  no  external  aids  such  as  LORANt,  this  study  concludes  that  successful  use  of 
TERrOM  depends  primarily  on  the  heading  accuracy  of  the  inertial  navigation  .system  cho.sen  for  the  aircraft . 

The  computer  programs  developcHl  for  this  study  can  he  used  with  larger  array  sizes  which  are  not  nece.s.sari ly 
.square  and  a longer  track  length,  permitting  the  .study  of  other  TP'.Rt'OM  configurations.  The  programs  which 
are  listed  in  appendices  A (Terrain  (lenerationl  and  H (P'light  Simulation'  provide  a llexihle  computational 
tool  for  the  determination  of  acceptable  accuracies  of  various  TERCOM  configurations. 
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APPENDIX  A 

TERRAIN  GENERATION  PROGRAM 


T-.r<.’(Xr.fMT,Ot'TruT,Tap-l  ) 

I ’ i*r[Or4  ^ ) 

I J=  1 

Fm  _ /\N-  J . 

L - li  • 

iAlL  <ANS>-T(£) 

SI  ■.NO  = ^'^. 

U'=  1 
N-'  X ■- 
XO=-.J. 

iip  10  A N'r<=i,to 
00  2 ')  I = 1 , I D 1 1 
"0  '■j  J-l.iOIM 
<I(  I . J)  -0. 
u3  CONIINU-’ 

IFi  TO.  'XJ.  r.)  r.O  TO  2U 

■{0  = 3. 

K0l=u. 

GO  TO 

O'.  K0  = :xp  (-Ix/TD 

1 = LX?(-S',3h  T (2.  ) / TO 
2I>  «CA  = X0/(1.“{01) 

Al,:)  =SijkT  (l.-2.*(('.'n)»*2»/(-{01*l.)) 

CAlL  '.ANr)?(KNl,SN2,SIGNO,F'(CAN,in) 

><G  i = '<  N 2 

00  TO  I=l,iniM,2 

OAlL  ^AN02<K■N1  ,f{N2,  SIGNO.FMtAN,  IC) 

A(  I , 1 ) i-IO’i^NTxSQ^T  ( 1 .-«0»*2)  •KNl 

A(  I ♦!  , 1)  =00»  A (I , T ) ♦SQ'^T  ( 1 . - “{O*  » 2 ) »RN2 

a;n  0 = A ( I ♦ 1 , 1 ) 

3 0 COixTlNUt 

DO  ->1  J = 2,I0IH 

call  >vANn'>lPNl,RN2,SH,NO,FM£AN,IO) 

Al  I , J)  =R0*A(  1,  J-1)  XSORT  ( 1 . - r<n»  » 2 ) ♦R  N1 
A(2,J):ACA*(A(1,J)*A(2,J-1>)+ACB*KN2 
DO  31  I = 3,ir'IM,2 

call  ■<ANU?(RNl,RN2.SIGNO,FMt:AN,  ID) 

A(  I , J)  = ACA»  ( A( 1-1 , J) .A ( I ,J-1) ) ♦ ACB*  ?N1 
A(  I *1  , J)  =ACA»  (A(I,JMa(I  + l.J-l))  ♦ACB»KN2 
31  CONTIN'Jf 

DO  3'5  I = 1,1  DIM 
00  3b  J=l,iniM 

A(l,J)=A(l,J)»(l./(l.* (FLOAT (J-l)/XO)»»NFX)) 
35  CONTINU' 

writ:  ( 1)  NCI  R,  ( ( A(  I ,J)  , 1=1 , IDIM)  ,J=  l.IDIM) 
loO  CONTINUl 

PRINT  102,NI'TR,l 

U2  format  (10X,»NCT£,  = »I2,5X,»>=»05. 7) 
call  RAN&tTIP) 

PRINT  103, P 

It?  FORMAT (13X,»P=»L15. 7) 

PtRIND  1 

STOP 

tNO 


I.^ 


iEL’Sp?  " 3''“”''  raAcnoBM 

OOPK  i OhMltitiAU  iV  iJpQ  -■ 


APPENDIX  B 

TERCOM  FLIGHT  SIMULATION  PROGRAM 


P-<Oo‘-’AM  TC-^  (INP'.n  , OUTPUT, T 4D£1  ) 
PINLNSION  4(&4,0'4)  ,S(b4)  ,43(0',f&4) 
01  If  N5I0N  UN  (64) 

•<0  = . ■T9')J23'9 
E=  . 7fr•.^67^ 

PNN=1. 09 
•NL  = 4? 

NOL  =b4 
NLS=NnL-NL 
9L5=FL0flT  (NL  5) 

3L  5 = 1 1 . 

CALL  94NSfT(F) 

PRINT  lO’.E 
PRINT  101 

ANUM  = S.NUN  = OMSO  = OHAn-0. 

00  7e  NO  = 1,10 
AOl =Sni=A00=S02=Q. 

KFAO(l)  NCTR,A 
PRINT  104, NC 
SUN10=0. 

SUN11:0. 

RA=20. 

PHI =0. 

PRINT  102, ^Hl 
PMl=(Pr)I/lflO.)*3.141tq26 
FMt  AN  = 0. 

IK= JK=64 
aN=FLOAT (IK* JK) 

SUN=0. 

BN2=FL04T(IK-1) 

00  10  1=1, IK 
00  10  J=1,JK 
SUN=SUNtA (I , J) 

10  CONTINUE 
Ave=SUN/BN 

suN=a. 

00  20  J*1,JK 
DO  20  1=1, IK 
SUN = SUM* (A(I,J)-AVE)**2 
20  CONTINUE 

SIGT=SQRT(SUM/BN) 

SI&NO=SIGT/R4 

10=1 

IF ( RO.Nt.O. ) GO  TO  22 
ROl=0. 

GO  TO  2S 

22  IF{R0.N£.1.)G0  TO  23 
R01=l. 

GO  TO  25 

23  SCN=-l./4LOG (RO) 

R01=iXP(-SaRT(2.) /SCN) 

25  ACA=R0/(1,*R01» 

AC8=SQRT( 1.-2.* ( (RO)** 2) /(R01*l . ) > 
call  RAN02(RN1,RN2,SIGNO,FNEAN, 10) 
RN3=RN2 
00  30  1=1,64,2 

call  RAN0?(RN1 ,RN2,S1GN0,FNEAN, 10) 
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THIS  PAOI  IS  BSST  QUALITY  rRACTlCABLJ 
TROH  COPY  FURtUSHkD  TO  WO 


■<N  t l*n  > ♦$  1«T  I 1 . •I^N^ 

SU’lt  0 

SU'<  1 IsKN  ( 1)  ♦RN  (!♦  1)  ♦SUMl  1 
AUT,U<A(I,1)  *ICN(  I> 

A3  ( 1 1 1 , 1 ) lA ( !♦ 1 , 1 1 tKN (I *11 
KN3  = -<N  (I  ♦ U 
o3  CONTINUE 

00  il  jsj,64 

call  «AN0r(ANl,‘f.»2,SlGN0,F«iLAN,  10) 

KtN(  1 ) = '^J»RN(  1)  tSO-lT  li  . -i«0»*2)  »«N1 
KN  ( 2)  iACA»(KN(  1)  »*^N(2)  ) < AC-)*i<N2 
SU'11  a = ’<Nl  1)  ••  >*<N(  2)  ••2*SU'110 
S0Mll  = .<Ml)*«Nl2)  ♦SUMll 
A0(  1,J)  :A(1,  J)  «i<N(  1) 

AO  t 2, J) =Al2, J) ♦KN (2) 

00  Jl  I<jtUHi2 

call  KAMr)2(KNl,KN2,SluN0,FH£AM,ID) 

^N(l)=ACA•(eN(l-^)♦^^(n)♦A':B•l<Nl  i 

KN  ( 1*1 ) s ACA*  (i<K(  I ) ♦kNl  I«l)  ) ♦ ACO  •KN2 
SU'H0=»NlI)»»2*r<NlI*l)*»2*SUMlJ 
SUMllsi<N(  1)  ♦KNlI  *1 ) ♦5LINI I 
AO  ( 1 , J)  >A  a , J)  ♦KN  ( I) 

AOH*i,j)iAlUl,J)t«N(I^l) 

Jl  CONTI  NO-: 

SUN:0. 

00  b2  J=1,JK 
00  12  I»1,1K 
SONiSUN^AOU  . J) 

*.2  COnIINUL 
SON! 

no  S'  jiifjK 

no  No  1 = 1,  IK 
A0( I , J) X A )(I , J)-S0N 

•'J  coNtiNun 

SJ:vNxS  )K1  iKUNlO/wOBf  . - ( SON  1 1 / h p oj, , )»»2) 

IF  ( SOt  VN.N,.,  J.  ) 1,0  to  4 0 

SNKxii.n  j 

CO  TO  Hi  I 

*.3  SNKxSIi.T/SOL  VN 

Hi  PKiNT 1 IJ.KOf VN.SlCT.SNR  , 

10  = 2 

00  2 T L X I , 1 C t 

call  KAN0r(«'Nl,KN2,SIliN0,FN£AN,  ID) 

TOx  (OLS'KND  ♦ 1. 

XOx  <hS.»KN2)  'IC. 

M 1 xro 

KJl xXO 
F N X ,1 . 

00  h2  Nxl.SL 
Xx**N»S,Nt)'Ml  ) ♦>0 
Va<>N»2JS(>Ml  ) Hio 

1 iHX  (T) 

J» IFIX (X) 

( »-►  lCA  T I I ) ) / ( x-Fl  OA  T IJ)  ) 

IF  l.)00  TO  hT 

o!  (A ( I,  j)-A  1 1 . jx  1 ) )•  (Float  ) jx  i ) -x) 


1 1 ■ 


K'  ) 


l.U' 


t .«•> 


t <•>} 


i-.r> 


ISO 


I 


>.'*  l A 1 1 M % J‘  U - AC  I , j«  U I *l  > -H.OAT  ( n ) 

Sell  •‘M  * ■.>.'•  A ( I , J»  1 ) 
wc'  T cl  ••4 

ai • 1 A c 1 ,,n-A  ( I ♦ 1 , ji ) •(FLOAT  c 1 ♦ I > -y  1 
(A(J»J,J«l)-A(J*l,J))»|X-FLaATljn 
Sell  «i> ( ♦ oa* A c I • 1 1 j I 
.1,  ('N»1'N»'’NN 
>•,'  CONTINJc 
SOM^O. 

110  '•> 

SUM •SUT*S  n ) 

*S  c'ONllNlir 
SUM  •bill/  NL 
00  se  :*l,NL 
Sll ) •S  11 1 -SUM 
*'  D c'O N T I N JL 

NL  .’•t'4-NL»l 
SliMK».U'Cl30i 
SlIMSO*  1.!  7 
00  hO  .(•liJK 
00  h.’  K»l,NL.' 

SIIM«0, 

SUMO*fl . 

SLIMNifl, 

00  I'.A  1»1,NL 
11  1*4-1 

SOMAiA  H 11,  J)  ♦SIIM4 
S.A  c,ONllNOL 

SOMHiSUMN/f  lOAT  (Nl 1 
00  hJ  1 l.Nl 
11* 1*4-1 

SUM.sliM*AaS<S(  D-AOlll  , J)*SUNS) 

SIIM.'»MIH.'*<S  ( n -Aa«  11  , J)  *S0N5)  • *^ 

1,1  I'ONTlNOf 

IF  I SOM.l.t  .SUMOGO  TOM 

SLIMS  •SUM 
4l  < *4 
4JM*J 

M IF  i sum:.i;i  .suMKO)  r.o  to  oo 

SUMS.’*  SON.* 

M <0*4 
4JS.’*J 
CONTlNUf. 

AOlST.SUSm  float  cKU-MK)  > ••?♦  (FLOAT  (KJl-KJRl  1 ••Tl 

S01Sl»SclST(  (FLOAT  (Kll-KiR.’t  » ••  .>  ♦ (FL  OAT  ( K Jl -K  JR?  I 1 ••?! 

OHINT  lOQ.AOIST.SOIST 

Ant •AniST*AtH 

SOI •soiST*sni 

AO.’«APIST»»,'*AO? 

SO.’  •SPIST*»?*SO? 

FO  I'ONIINUF 
AO  I • AO  1/ 10. 


SOI *501/10. 

Ati.’^An.’/lc). 

sic*sn.'/iP. 

F'SINT  l.’c),A01,AO?,Sin,SO.' 
7.'  l'ONTINU! 


.jcS* 
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REMIND 

PRINT 

PRINT 

PRINT 

call  R 
PRINT 
STOP 

100  FORMAT 
ICl  FORMAT 

102  FORMAT 

103  FORMAT 
ION  FORMAT 

105  FORMAT 

106  FORMAT 
lOF  FORMAT 
108  FORMAT 
lOR  FORMAT 

110  FORMAT 

111  FORMAT 

112  FORMAT 

113  FORMAT 
1»S/N»» 

IIM  FORMAT 
120  FORMAT 
lElO.R) 
END 


103, RA 
112, RO 
11<«,NI 
ANOET (00) 

107,00 

(16X,F6.2,27X,F6.2) 

(13X,»MAO  ERROR* ,24X,»MS0  tRROK*) 

(1X,*PMI»*,F4.1) 

(IX, ’SIGMA  T / SIGMA  N s»,F4.1) 

(1X,»NCTR=*, 13) 

(lX,*CORKtCT  1. 0. AT*, 213,  lOX, 213, IX, •MISTAKEN  FOR*, 213) 

(IX, ’CORRECT  I.O. AT*, 213, lOX, ’CORRECT  I. 0. AT*, 213) 

(IX, ’RANDOM  NUMBER=*,E15. 7) 

(IX, 213, IX, ’MISTAKEN  F OR* , 21 3 , 5X , ’CORRECT  I.O. AT*, 213) 

(IX, 213, IX, ’MISTAKEN  FOR* , 21 3 , 5X , 2 1 3 , 1 X , * MI  ST AKFN  FOR*, 213) 
(7X,*NA0*,32X,*MS0*) 

(lX,*AVt.  MAO  ERRORS*, IX, F5,l,llX,*A«t.  MSO  ERROR”  , 1 X , F5 . 1 ) 
(1X,*R0**,F5.2) 

(bOX, ’NOISE  SIGMA**, E15. 7, IX, ’SIGNAL  S I GM A* * , E 1 5 . 7 , 

,F5.2) 

(IX, ’TRACK  LENGTH  **,I3) 

(1X,*A01=*,£10.4,1X,*A02**,E10.4,IOX,*S01»*,E10.4,1X,*S32**, 


IS  ?AQE  IS  BSST  QUALITY  PRAOXI^ABIJI 
OOfl  fURBlSHEi)  IV  DD.C 


K4Nl\’(K'Nlf«Ni,SI&N0t  FHt*N*IP) 

y«-?4Nt  invi**) 

If  tin. fa. 1100  10  •> 

RNJil 

uO  TO  to 

cc«sa«n-J.»»oo&(x>  »*siGNO 

><N  I ■ CC 'C OS  ( f 3 1 J • V ) ♦ F *»F *N 

i<N;»CC»SIN(t.,  J«3l3»f  I yf  Hf4N 
KSTUKN 
tNO 
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